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We have reexamined the problem of n — n oscillations for ultra-cold neutrons (UCN) confined 
within a trap. It is shown that for up to 10 3 collisions with the walls the process can be described 
in terms of wave packets. The exact equation is derived describing the time-evolution of the n- 
component for an arbitrary number of collisions with the trap walls. The n-component grows 
linearly with time with enhancement factor depending on the reflection properties of the walls. 
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For quite a long time physics beyond the Standart Model continues to be an intriguing subject. Several reactions 
Q, which may serve as signatures for the new physics have been discussed. One of the most elegant proposal is to look 
for n — n oscillations (see also Q). There are three possible experimental settings aimed at the observation of 
C/0 ! this process. The first one is to establish a limit on nuclear instability since n produced inside a nucleus will blow it 
H ' U P- The second one is to use a neutron beam from a reactor. This beam propagates a long distance to the target in 
which possible n component would annihilate and thus be detected. The third option which we discuss in the present 
paper is to use ultra-cold neutrons (UCN) confined in a trap. The main question is to what extent the generation 
of the n component is reduced by the interaction with the trap walls. This subject was addressed by several authors 
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... ,. ..„ „. ,. j. A thorough investigation of the problem is in our opinion still lacking. 
| First of all the clear formulation of the problem of n — n oscillations in a cavity has been hither to missing. Two 
different approaches were employed without presenting sound arguments in favor of their applicability and without 
\Q 1 tracing connections between them. 

In the first approach 0, 0| n — n oscillations are considered in the basis of the discrete eigenstates of the trap 
potential taking into account the splitting between n and n levels and n annihilation. The density of the trap 
eigenstates which is proportional to the macroscopic trap volume is huge and the states cluster together extremely 
thickly. However these arguments are not enough to discard the discrete states approach since n — n mixing parameter 
Qh| is much smaller than the distance between adjacent levels — see below. The true reason due to which the discrete 
JL . eigenstates approach is of little physical relevance is the following. The spectrum of the neutrons provided to the trap 
qj' by the source is continuous and certain time is needed for the rearrangement of the initial wave function into standing 
waves corresponding to the trap eigenstates. As will be shown below this time interval is of the order of /3-decay time 
so that the standing waves regime being interesting by itself can hardly be reached in the real physical situation. 

The second approach 0,1a, uj treats the neutrons and antineutrons inside a trap as freely moving particles which 
undergo reflections from the trap walls. Collisions with the walls result in a reduction of the n component compared 
to the case of the free-space evolution. This suppression is due to two factors. The first one is the annihilation inside 
the walls. The second one is the phase decoherence of n and n components induced by the difference of the walls 
potentials acting on n and n. Reflections of antineutrons from the trap walls were for the first time considered in 
0]. The purpose of that paper was to investigate the principal possibility to observe n — n oscillations in a trap and 
the authors estimated reflection coefficient for antineutrons without paying attention to the decoherence phenomena. 
Only single collision with the trap wall was considered in |3) . A comprehensive study of n — n oscillations in a trap 
was presented in Decoherence and multiple reflections were included as well as the influence of the gravitational 

and magnetic fields. The approximate equation for the annihilation probability after N collisions was obtained in Q 
(equation (3.8)) which coincides with the exact formula (59) of the present paper when N ^> 1. As we show below 
(see Fig. nj the A^-independent asymptotic regime settles at N > 10. 

The derivation of the exact equation for the annihilation probability for arbitrary number of collisions is in no way 
the only purpose of the present work. We already mentioned the problem of the relation between the eigenvalue 
and the wave packet approaches. Within the wave packet approach some basic notions such as the time between 
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successive collisions and the collision time itself can be defined in a clear and rigorous way. Another question within 
the wave-packet formalism is the independence of the reflection coefficient on the width of the wave packet and the 
applicability of the stationary formalism to calculate reflections from the trap walls. These and some other principal 
points are considered for the first time in detail in the present paper. 

Let us also mention that an alternative approach to the evaluation of the reflection coefficients for n and n was 
outlined in Q . Calculations presented there are based on the time-dependent Hamilton formalism for the interaction 
of n and n with the trap walls. This subject remains outside the scope of the present paper. 

The paper is organized as follows. In Section II we remind the basic equations describing n — n oscillations in free 
space. Section III is devoted to the optical potential approach to the interaction of n and n with the trap walls. In 
Section IV we analyze the two formalisms proposed to treat n — n oscillations in the cavity, namely box eigenstates 
and wave packets. In Section V the reflection from the trap walls is considered. Section VI contains the main result 
of this work which is the time dependence of the n component production probability. In Section VII conclusions are 
formulated and problems to be solved outlined. 

II. OSCILLATIONS IN FREE SPACE 

We start by reminding the basic equations describing n — n oscillations in free space. The phenomenological 
Hamiltonian is a 2 x 2 matrix in the basis of the two-component n — n wave function (we set % = 1) 



H jl = {H j -i^-)8 jl +e{a x ) jl , (1) 



where j, I — n,n; Hj = k 2 /2m — ujB, fj,j is the magnetic moment, B is the external (e.g. the Earth) magnetic field, 
Tp is the /3-decay width, e is the n — n, mixing parameter (see below), cr x is the Pauli matrix. Assuming the n and fi 
wave functions to be plane waves we write the two-component wave function of the n — n system as 

= (Hf) e**. (2) 



Then the evolution of the time-dependent part of ^(Xji) is described by the equation 



The difference between E n and E n due to the Earth magnetic field is 



(3) 



lu = E n - E n = 2\n n \B w 6 x 1CT 12 eV. (4) 
Diagonalizing the matrix in (3) we find tp„(t) and ipn(t) m terms of their values at t = 0: 

i>n(t) — (Vn(0)(cos^ + ^sin^) -^(0)-sinz/t) exp[--(iQ + T p )t] (5) 

ipn(t)= f-^ n (O)^sinz4 + <0 n (O)(cosz4 - ^sin^)^ exp[~(iSl + Tp)t], (6) 
where = E n + E n , v = (uj 2 /4 + e 2 ) 1 / 2 , w = E„ - E n . In particular if V'n(O) = 1, ipn(0) = 0, one has 

I^WI 2 = JT^2 exp(~r^)sin 2 (iV^+4^t). (7) 
The use of this equation to test fundamental symmetries is discussed in . 
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Without magnetic field, i.e. for u> = 0, and for i«f 1 equation (7) yields 

|Vfi(*)| 2 «e 2 t 2 exp(-iy). (8) 

This law (for t <C T^ 1 ) has been employed to establish the lower limit on the oscillation time r = e _1 . According to 
the ILL-Grenoble experiment ^3] 

t > 0.86 x 10 8 s. (9) 

The corresponding value of the mixing parameter is e w 10~ 23 eV. This number will be used in obtaining numerical 
results presented below. 

The Earth magnetic field leads to a strong suppression of n — n oscillations. With the value of w given by (4) 
equation (7) leads to the following result: 

4e 2 

Wn{t)\ 2 « -j exp(-r>i) sin 2 t/rn « 10~ 23 sin 2 i/r B , (10) 

where tb = (!/■*« I-B) -1 «2x 10~ 4 s. In what follows we shall assume that magnetic field is screened. 
For w = but for arbitrary initial conditions equations (5-6) take the form 

Mt) = (^„(0)coset-^(0)sme*)exp[-(lE+-^)t] (11) 
^n(t) = (-iV'n(0)sinet + V»i(0)coset)exp[-(i J E;+ -£)t], (12) 
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where £7 = E n = E n . 



III. OPTICAL POTENTIAL MODEL FOR THE TRAP WALL 

We remind that neutrons with e nerg y E < 10 _7 eV are called ultra-cold. An excellent review of UCN physics was 
given by I. M. Frank [IjJ (see also [jj). 

A useful relation connecting the neutron velocity v in cm/s and E in eV reads 

t)(cm/s) = 10 2 [10 9 £(eV)/5.22] 1/2 . (13) 

For E = 10~ 7 eV the velocity is v ~ 4.4 x 10 2 cm/s. 

A less formal definition of UCN involves a notion of the real part of the optical potential corresponding to the trap 
material (see below). Neutrons with energies less than the height of this potential are called ultra-cold. The two 
definitions are essentially equivalent since as we shall see for most materials the real part of the optical potential is of 
the order of 10 ~ 7 eV. 

Our main interest concerns strongly absorptive interaction of the h component with the trap walls. Therefore wery 
weak absorption of UCN on the walls (Til Il2 | will be ignored. Due to complete reflection from the trap walls UCN 
can be stored for about 10 3 s (/3-decay time) as was first pointed out by Ya. B. Zeldovich 13] . 

To be concrete we consider UCN with E — 0.8 x 10 -7 eV which corresponds to v = 3.9 x 10 2 cm/s (see (13)), 
k = 12.3 eV and de Broglie wave length A « 10~ 5 cm. In the next Section we shall describe UCN in terms of the 
wave packets and hence the above values should be attributed to the center of the packet. 

Interaction of n and n with the trap walls will be treated in terms of energy-independent optical potential. The 
validity of this approach to UCN has been justified in a number of papers, see e.g. 0, 0, Q|. There is still an 
open question concerning the discrepancy between theoretical prediction and experimental data on UCN absorption. 
Interesting by itself this problem is outside the scope of our work since as already mentioned absorption of neutrons 
may be ignored in n — n oscillation process. The low-energy optical potential reads 

U jA = —Na jA , (14) 
m 
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where j = n, n; m is the neutron mass, N is the number of nuclei in a unit volume, clja is the j — A scattering length 
which is real for n and complex for n. For neutrons the scattering lengths a n A are accurately known for various 
materials 0|. For antineutrons the situation is different. Experimental data on low-energy n — A interaction are 
absent. Only some indirect information may be gained from level shifts in antiprotonic atoms. Therefore the values 
of dfiA used in |3, |6J, |& LL5J as an input in the problem of n — n oscillations are similar but not the same. We consider 
as most reliable the set of anA calculated in ^| within the framework of internuclear cascade model. Even this 
particular model leads to several solutions. Therefore the one we have chosen for 12 C (graphite and diamond) may 
be called "motivated" by Ref. dfj. To get a feeling on the dependence upon the material of the walls and to compare 
our results with that of 0] we also made calculations for Cu. Scattering lengths for Cu are not presented in ^| and 
we used a solution proposed in Q ■ Thus our calculations were performed with the following set of n — A scattering 
lengths: 

a»ic = (3 - il) fm, a nCu = (5 - i0.5) fm. (15) 

The scattering lengths for neutrons are [l2T ]: 

a nC = 6.65 fm, a nCu = 7.6 fm. (16) 

The concentrations of atoms N entering into (14) are: iVc(graphitc) = 1-13 x 10~ 16 far 3 , Afc (diamond) = 1-63 x 
lCT 16 fm~ 3 , iV Cu = 0.84 x 1CT 16 far 3 . Then according to (14) the optical potentials read 

E4c( g r) = 1.95 x 1CT 7 eV, 
^c(diam) = 2.8 x 1(T 7 eV, (17) 
U nCu = 1-66 x 1CT 7 eV; 

C4c( g r) = (0.9 - i0.3) x 1(T 7 eV, 
t4c(diam) - (1.3 - i0.4) x 10- 7 eV, (18) 
U nC n = (2 - i0.2) x 10" 7 eV. 

In the present work we consider particles (n and n) with energies below the potential barrier formed by the real part 
of the potential. For n and 12 C the limiting velocity is v = 4.15 x 10 2 cm/s. 

IV. WAVE PACKET VERSUS STANDING WAVES 

It is convenient to use for the optical potentials (17) and (18) the shorthand notation 

U^Vj-iWjSjn, (19) 

where j = n,n while the wall material is not indicated explicitly. We consider the following model for the trap in 
which n — fi oscillations may be possibly observed. Imagine the two walls of the type (19) separated by a distance 
L ~ 10 2 cm, i.e. the one-dimensional potential well of the form 

Uj{x) = {6{-x — L) + 0(x)} {Vj - iW 3 5 jn } , (20) 

with 9{x) being a step function. Our goal is to follow the time evolution of the n component in such a trap assuming 
that the initial state is a pure n one. 

The first question to be answered is how to describe the wave function of the system. Two different approaches 
seem to be feasible and both were discussed in the literature 0.0. Hi- The first one is to consider oscillations occuring 
in the wave packet and to investigate to what extent reflections from the walls distort the picture as compared to the 
free-space regime. The second one is to consider eigenvalue problem in the potential well (20), to find energy levels 
for n and n, and to consider oscillations in this basis. Due to different interaction with the walls the levels of n and 
n are splitted and n levels acquire annihilation widths. 
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At first glance this approach might seem inadequate since in a trap with L ~ 10 2 cm the density of states is very 
high, the characteristic quantum number corresponding to the UCN energy is very large, the splitting SE between 
adjacent n-levels (or between the levels of n and n spectra) is extremely small. The values of all these quantities will 
be given below and we shall see that SE < 10~ 14 eV. However, this approach can not be discarded without further 
analysis since the n — n mixing parameter e s=s 10~ 23 eV is much smaller than SE. 

To understand the relation between the two approaches note that the initial conditions correspond to a beam of 
UCN provided by a source. The momentum spectrum of UCN depends on the concrete experimental conditions. In 
order to stay on general grounds and at the same time to simplify the problem we shall assume that the UCN beam 
entering the trap has a form of a Gaussian wave packet. Suppose that at t = the center of the wave packet is at 
x = x so that 

VfcOr, t = 0) = (Tra^-V^exp (- {X ~J o) + ikx^j , (21) 

where a is the width of the wave packet in coordinate space. The normalization of the wave function (21) corresponds 
to one particle over the whole one-dimensional space: 

+ oo 

|2 



J dx\tp k (x,i 



,t = 0)\ z = l. (22) 

For E = 0.8 x 1CT 7 eV and the beam resolution AE/E = 10 3 one has 

k = 12.3 eV, a = 3.2 x 10~ 3 cm. (23) 
The width of the wave packet (21) increases with time according to 



a = a 



t 



1/2 



t/ma (24) 



and for t ~ 10 3 s becomes comparable with the trap size L. In order that the wave hitting the wall and the reflected 
one to be clearly resolved the condititon a' /v <C tl, or a' <C L has to be satisfied, where tl <~ 1 s is the time between 
two consecutive collisions with the trap walls. Reflection of the wave packet from the walls is considered in detail in 
the next Section. Here we will show that t ~ 10 3 s is the characteristic time needed for the rearrangement of the 
initial wave packet into the stationary states of the trapping box. 

Consider the eigenvalue problem for the potential well (20). Parameters of the potential (20) for neutrons are 
V n ~ 2 x 10~ 7 eV, L w 10 2 cm. The number of levels is 

M « „ 10 8 /7r . ( 25 ) 

According to (23) the center of the wave packet (21) has a momentum k = 12.3 eV which corresponds to a state 
with a number of nodes j~2x 10 7 and kjL w 6 x 10 7 3> 1. Positions of such highly excited levels in a finite depth 
potential are indistinguishable from the spectrum in a potential box with infinite walls. Thus 

(fij(x) rts J —sinujjX, ujj = -j-. (26) 

The wave functions (26) describe semi-classical states with j ' ^> 1 in a potential well with sharp edges. The "frequency" 
u)j is very high as compared to the width of the wave packet in momentum space, u 3 ; ~ 6 x 10 5 cm -1 ^> v = l/(\/2a) ~ 
2 x 10 2 cm -1 . This means that the wave packet spans over a large number of levels. To determine this number note 
that the distance between adjacent levels around the center of the wave packet is SE = Ej+i — Ej <~ 10~ 14 eV. 
The highly excited levels within the energy band AE = 10~ 3 E ~ 10~ 10 eV corresponding to the wave packet (21) 
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are to a high accuracy equidistant as it should be in a semi-classical regime. The number of states within AE is 
Aj = AE/SE ~ 10 4 and their density in momentum space is 



p(io) = aAj c± L/tt ^ 10 6 cV 1 . (27) 

Now we can answer the question formulated at the beginning of this Section, namely whether one should describe 
n — n oscillations in the trap in terms of the wave packet or in terms of the stationary eigenfunctions. At t = the 
wave function has the form of the wave packet (21) provided by the UCN source. Due to collisions with the trap walls 
transitions from the initial state (21) into discrete (or quasi-discrete for n) eigenstates (26) take place. 

The time-evolution of the initial wave function (21) proceeds according to 

4>(x,t) = J dx'G(x,t;x',0)ip k {x',0), (28) 

where G(x,t;x',0) is the time-dependent Green's function for the potential well (20). Making use of the spectral 
representation for G, one can write 

ip{x,t)=^2e- iE ' t ipj(x) f dx'<p*(x')Mx',0). (29) 

3 

In semi-classical approximation the distance between the adjacent levels is SE — tt/tl- Therefore one may think 
that at t = tl, i.e. already at the first collision the neighboring terms in (29) would cancel each other. However this 
is not the case. Indeed, 



e -iEit 



<p j+1 (x)e- iE * +lt + ipj{x)e~ lE ^ = —== e iuj * x {l + e 4 f (^-«*)) _ e -*"i*(i + e -*£(*-H'*)) 

iv2L L 



Therefore at x = ±vt there is a constructive interference either in the first or in the second terms correspondingly. 
This holds true with the account of the whole sum of terms in (29) and hence one can pass from summation to 
integration in (29) . The wave function overlap entering into (29) can be easily evaluated provided the center of the 
wave packet xq is not within the bandwidth distance a' from the trap walls. The overlap is given by the following 
integral: 

o 

f dx>v*(x')Mx',0)^ (2 ^ o) i/ 2 / dx' Cxp [- {x ' 2 J o)2 + l (k-co J )x' 

^ ^ (3 °) 

i [ , „ ( (x"+x Q ) 2 „\ 



o 



At this step we have omitted the exponent with high frequency (k + u)j). Next take (x" + xo)/(\/2a) as a new variable 
and resort to the assumption \xo\ ^> a, L — \xq\ 3> a (we remind that xq is negative since — L < x < 0). The result 
reads 

J dxV*(x')^(^,0)^*(^) 1/2 exp^-y(A : -^) 2 +z(A : -^)xo^ . (31) 

Corrections to (31) are of order of a/L. Now comes frequency summation in (29). This summation can be substituted 
by integration over u> since the density of semi-classical states p(u>) is very high. In this way we arrive at 

, , . 1 ( , aix.i) . , a(—x,t) ,\ , , 



6 



(x - in) 2 k 2 a 2 t 2 

a(x, t) = (x-x - v t) 2 - it- ^ 2ik a 2 (x - x ) + i— — t + 2ik x (a 2 H 5-A (33) 

ma m m' : a z 

The second term in (32) describes the reflected wave packet — see the next Section. According to (21), (28) and (32) 
all that happens to the wave packet in the trap is broadening and reflections. This is true at least during some initial 
period of its life history. How long does this period last? The answer to this question may be obtained estimating 
the accuracy of performing frequency integration instead of summation over discrete states in (29) . 

To estimate the time scale for the rearrangement of the initial wave packet (21) into the trap standing waves (26), 
it is convenient to introduce the difference 5ip(x,t) = "4> S um(x,t) — ipi n t(x,t) between the "exact" wave function (29) 
and the approximate integral representation (32). As soon as 

Sw(t) = J dx (|t/w| 2 - \^nt\ 2 ) = 2 J dx^ m M) « 1, (34) 
one can consider oscillations as proceeding in the wave packet basis. The following estimate holds 



6ip(x, t)=^2 fi^n) - / dup(uj)f(uj) = / dujp(uj)(f(uj) - f(u n )) 

- - ^ J rfwp(w)/'(w n )(w-w n ) = -^/'WK+i -<*>„); 



here 



From (36) one gets 

g(u>) = i(x — xo — vt) — (fco — uo)a 2 . 
Since f(u) is a narrow Gaussian peak, we can substiute g{u) by g{k ), and then (35) results in 



From (34) and (38), we get 



+00 

t t 



maL 10 3 s ' 



(35) 



U/na f a 2 2 u! 2 ■ 

= " ~2I? ° XP V — 2^ * - i—t + ilj(x - x ) + ikxo ) . (36) 



(37) 



7T 

5ip(x,t) ~ —{x - x - v t)il>i n t(x,t). (38) 



(39) 



where a 1 is given by (23). 

Roughly speaking, the time t <~ 10 3 s needed for the neutron wave function to rearrange into the trap eigenstate 
is comparable to the neutron life-time and neutron would rather "die" than adjust to the new boundary conditions. 
Therefore in what follows the wave packet formalism will be used. Some additional subtleties arising from the 
quantization of the levels in the trapping box will be discussed in Section VII. 

V. REFLECTION FROM THE TRAP WALLS 

Consider again a one-dimensional trap (20). Let the particle moving from x = —00 enter the trap at t = through 
the window at x = — L. At t = tl it will reach the wall at x = 0, the n component will be reflected from the wall 
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while the n component will be partly reflected and partly absorbed. The wave packet describing the interaction with 
the wall has the form 



+00 

ip{x,t) = 7r_3/4 \/ r f /* dki> j (k,x)e-xp(~(k-k ) 2 + tL{k - k ) ~ i^k 2 ) , (40) 



where j = n,n and 



il>j(k,x) = e lkx + R{k)e- lkx = e' kx + Pj (k)e^^e- lkx . (41) 



For the n component p„[k) = 1 since we neglect very weak absorption of neutrons at the surface. The integral (40) 
with the first term of (41) is trivial. In order to integrate the second term of (41) note that due to the Gaussian 
form-factor with ako ~ 10 3 3> 1 the dominant contribution to the integral (40) comes from the narrow interval of k 
around ko- Expanding Rj(k) at k — k and keeping the leading term we get 

Rj(k) ~P3(ko)e^ {k °\l + iti(k )(k~ko) + 6^?^(k-ko)} , 

3 Pj(k ) (42) 

The validity of the last step for n will become clear from the explicit expression for pn(k) and cpn{k) presented below. 
Now integration in (40) can be easily performed with the result [TH ] 



a re ji(x, t) 



2 \ J I ' 



(43) 



(x + L) 2 k 2 a 2 t 2 
a inc (x, t) = (x + L- v t) 2 - i V _J - 2ik a 2 {x + L) + i-^-t + 2ik L(a 2 H — ), (44) 



a re fi(x, t) = a mc (-x + (j)' , t) + 2ik 4>'{a 2 H — ). (45) 

From (43)-(45) we see that the essence of R(k) in the wave packet formalism is the same as in the time independent 
approach. Therefore imposing standard boundary conditions at x = we get the reflection coefficients 

B J (k)=p J (k)e^ = ^^-, (46) 

K n = [2m{V n -E)] l '\ 
n fl = [2m(Vn - iW fl - E)] 1 ' 2 = 4 - in'i, 

Pn = ' pfi = (fcW+W (49) 

In particular for 12 C (graphite) 

p = 0.56, 9 = cj> n - <j> n = 0.72. (50) 
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The first term in the right-hand side of (45) may be written as \—x + L — Vo(t — 4>' /vq)] 2 . Hence the collision time 
or time-delay is 0, 0] 



2 JYI 

Tj.eoll = <f>'j(ko)/v = Re— . (51) 

For neutrons, i.e. for real n n , equation (51) gives the well-known result T n co n = [E(V n - E)]- 1 ' 2 . This result is m 
line with the naive estimate r njCO u ~ 1/vq ~ 10~ 8 s Q, where I < A is the penetration depth. 
For 12 C (graphite) equation (51) yields 



r n ,coii = 0.7 x 10~ 8 s, 7H,coU = 1-1 x 1(T 8 s. (52) 

Equations (43)-(45) supplemented by the above inequality allow to follow the time evolution of the beam inside 
the trap. Imagine an observer placed at a bandwidth distance from the wall, i.e. at X = —a. According to (43)— (45) 
such an observer will conclude that at times t < tl — r a the incident wave (the first term in (43)) dominates, while 
at t > tl + r a the reflected wave prevails. With this splitting of the time interval around Ntl, N = 1, 2, . . . in mind 
we shall use the notations (Nt^ — ) and (Ntl+) for the moments before and after the iV th collision. Thus, we can 
calculate n production rate since we have a rigorous definitions of the collision time and the time interval between 
the two subsequent collisions. 



VI. ANNIHILATION RATE IN A TRAP 



Now we can inquire into the problem of the time-dependence of n production probability. In free space it is given 
by 1^(4) | 2 = e 2 t 2 (see (2)) while in a trap with complete annihilation or total loss of coherence at each collision it 
has a linear time dependence |V'«(*)| 2 — e 2 TE,i 01< 

To avoid cumbersome equations and in view of the fact that we consider time interval t <C T^ 1 we omit exp(— Tpt) 
factors. Production of n during the collision can be also neglected 0. Difference in collision times (52) for n and n 
may be ignored as well. In the previous Section we have seen that the interaction of the wave packet with the wall is 
described in terms of the reflection coefficients (46). (2fJ| 

Assume that at t = a pure n beam enters the trap at x = —L. After crossing the trap, i.e. at t = (tl—), the 
time-dependent parts of the wave functions are given by (12) [we state this because although the Gaussian form factor 
in (43) also depends on time, the corresponding terms in the time-dependent Shrodinger equation are of order of 
1/afeo as compared to the derivative of the exponent exp(— iEt); also note that the form factors are the same for n 
and n up to constant multiplier]: 

ipn{T L -) = cos(er L ) exp(-iET L ), 
Mtl-) = sin(er L ) eaq?[-i(Er L + tt/2)]. 

After the first reflection at t = (tl+) we get 

il>n{TL+) = cos(eTi) exp[—i(ErL - 4> n )], 
ipn{TL+) = Pa sin(c7x) exp[-i(Er L - <pn + tt/2)]. 

Evolution from t = (t£+) to t = (2r^— ) again proceeds according to (12) 

4>n = I sin(2er L ) (l + pe l9 ) exp[-i(2Sr L - fa + tt/2)] 

* (55) 

~ er L (1 + pe ld ) ex V [-i{2Er L - fa + tt/2)], 
where 6 — fa — fa is the decoherence phase and p = p fi . Now the answer for iP(Ntl-) seems evident: 

1 - o N e me 

*fo{Nr L -) - er L ig exp[-i(NET L -fa+ tt/2)}. (56) 

1 — pe la 
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This conjecture is easy to verify by mathematical induction. For t = (2tl — ) the result has been derived explicitly — 
see (55). Evolving (56) through one reflection at t = Ntl and free propagation from t — (Ntl+) to t = [(N + 1)tl— ] 
we arrive at (56) with (N + 1) instead of N. This completes the proof. 
Thus the admixture of n before the N th collision, i.e. at t — Nt^— is 



I, , w v,a 2 2 l + P 2jV -2p N cas NO 

\MNr L -)\ =er L 1 + p2 _ 2pcosd (57) 
The annihilation probability at the j th collision is 



Pa(j) = (l-p 2 )\MjTL-)\ 2 . (58) 

Hence the total annihilation probability after N collisions is 



JV 



P a (N) = (l-p 2 )J2\MkTL)f 



k=l 

„2/-i „2N\ „„„Q n „N\ 



e 2 T 2 (l~p 2 ) ( p 2 (l-p 2N ) cos9-p-p N [cos(N+l)O + pcosN0] 



N + ^ " 2 P 



1 + p 2 - 2pcos9 \ 1-p 2 l + p 2 -2pcos9 



(59) 



After several collisions the terms proportional to p N , p 2N and p N+1 may be dropped since p ~ 0.5 (see (50)). Then 
(59) takes the form 



PJN)~ fi^ (n(1~ P 2 ) + 1 ^ V (60) 

aK ' l + p 2 -2pcos6 \ y F! l + p 2 -2pcos6J K ' 

Three different regimes may be inferred from (60). For very strong annihilation, i.e. p -C 1 



P a (N) = e 2 r 2 N = e 2 T L t. (61) 
For complete decoherence at each collision, i.e. for 8 = ir 



2 2 I A A - P , P( 2 - P)\ ^, 1 - P 2 



P a (N) = e 2 r 2 L ( + ) ^ ^-e 2 r L t. (62) 



For the (unrealistic) situation when 9 = 



For the values of p and 9 corresponding to the optical potentials (17)-(18) the quantity Q a {N) = (e 2 r|iV) _1 = 
(e 2 Ti J t)~ 1 P a (N) calculated according the exact equation (59) is displayed in Fig.Q] This figure shows that the linear 
time dependence settles after about 10 collisions with the trap walls. The asymptotic value of Q a (N) which may be 
called the enhancement factor is 1.5 ~ 2 depending on the wall material. 

Proposals have been discussed in the literature |6J, yj| to compensate the decoherence phase 9 by applying the 
external magnetic field. Suppose that in such a way the regime 9 = may be achieved. Also suppose that one can 
vary the reflection coefficient p in a whole range by varying the trap material. For such an ideal situation we plot in 
Fig. [2 the quantity N e ff(p) defined as 



P a {N)=e 2 r 2 N eff {p). (64) 

Thus defined N e ff(p) obviously depends also on the number of collisions N and in Fig. |2the results for N = 10 
and N = 50 are presented. This figure shows what can be expected from the trap experiments in the most favorable 
though hardly realistic scenario. 
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VII. CONCLUDING REMARKS 



We have reexamined the problem of n— n oscillations for UCN in a trap. Our aim was to present a clear formulation of 
the problem, to calculate the amplitude of the n component for arbitrary observation time and for any given reflection 
properties of the trap walls. We have shown that for physically relevant observation time (i.e. for the time interval 
less than /3-decay time) the process of n — n oscillations is described in terms of the wave packets while the standing 
waves regime may settle only at later times. By calculating the difference between n and n collision times the new 
light has been shed on the decoherence phenomena. For the first time an exact equation has been derived for the 
annihilation probability for an arbitrary number of collisions with the trap walls. In line with the conclusions of the 
previous authors on the subject this probability grows linearly with time. We have calculated the enhancement factor 
entering into this linear time dependence and found this factor to be 1.5 2 depending on the reflection properties of 
the wall material. 

Despite the extensive investigations reviewed in this article and the results of the present paper the list of problems 
for further work is large. The central and most difficult task is to obtain reliable parameters of the optical potential 
for antineutrons. The beam of n with energy in the range of 10" 7 eV will be hardly accessible in the near future. 
Therefore work has to be continued along the two lines mentioned above — to deduce the parameters of the optical 
potential from the level shifts in antiprotonic atoms and to construct reliable optical models which can be confronted 
with the available experimental data on n-nuclear interaction at higher energies. In the forthcoming publication we 
plan to present numerical calculation of the time evolution of the wave packet into standing waves as well as to discuss 
some features of n — n oscillations in the eigenfunctions basis which were not discussed in Ref . [5| . Another task is 
to perform calculation for the concrete geometry of the trap an realistic spectrum of the neutron beam. To this end 
one needs an input corresponding to a concrete experimental setting. 

The authors are grateful to Yu. A. Kamyshkov for interest to the work and discussions. Useful remarks and 
suggestions from L. N. Bogdanova, F. S. Dzheparov, A. I. Frank, V. D. Mur, V. S. Popov, G. V. Danilian are 
gratefully acknowledged. 

One of the authors (B. O. K.) expresses his gratitude for financial support to V. A. Novikov, L. B. Okun, INTAS 
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Figure 1: Plot of Q a (N) = (e 2 T L t)~ 1 P a (N) dependence versus N. Solid line corresponds to 12 C (graphite), dashed one - to 
12 C (diamond), dotted one - to Cu. 
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Figure 2: Plot of N e ff versus p dependence. Solid line is for the number of collisions N = 50, dashed one corresponds to 
N = 10. 
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